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ABSTRACT 

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star for- 
mation rates in galaxies is investigated for a large sample of galaxies observed with 
the SPIRE and PACS instruments on Herschel as part of the HerMES project. We 
build flux- limited 250 /iin samples of sources at redshift z < 1, cross-matched with 
the Spitzer /MIPS and GALEX catalogues. About 60 % of the Herschel sources are 
detected in UV. The total IR luminosities, LrR, of the sources are estimated using a 
SED-fitting code that fits to fluxes between 24 and 500 /xm. Dust attenuation is dis- 
cussed on the basis of commonly-used diagnostics: the Ltr,/Ltjv ratio and the slope, 
/3, of the UV continuum. A mean dust attenuation Auy of ~ 3 mag is measured in 
the samples. Lir./Ltjv is found to correlate with Ljr. Galaxies with Ljr > lO 11 !.© 
and 0.5 < z < 1 exhibit a mean dust attenuation A\jy about 0.7 mag lower than 
that found for their local counterparts, although with a large dispersion. Our galaxy 
samples span a large range of /3 and Lir/Ltjv values which, for the most part, are 
distributed between the ranges defined by the relations found locally for starburst and 
normal star-forming galaxies. As a consequence the recipe commonly applied to local 
starbursts is found to overestimate the dust attenuation correction in our galaxy sam- 
ple by a factor ^2 — 3. The SFRs deduced from Ltr are found to account for about 
90% of the total SFR; this percentage drops to 71% for galaxies with SFR < lM yr _1 
(or Lm < 1O 1O L0). For these faint objects, one needs to combine UV and IR emissions 
to obtain an accurate measure of the SFR. 

Key words: galaxies: evolution-galaxies: stellar content-infrared: galaxies-ultraviolet: 
galaxies 



1 INTRODUCTION 

Far-infrared (IR) and ultraviolet (UV) luminosities are com- 
monly used to estimate the current star formation rate 
(SFR) in galaxies since both emissions are expected to come 
from young stars. By combining observed IR and UV lu- 
minosities one can therefore make an energetic budget and 
derive an accurate measure of the total SFR in star-forming 
galaxies (e.g. Iglesias-Paramo et al.|2004 Elbaz et al.| 2007). 
It is often the case that only UV or IR data exist, however, 
and so the question remains of how reliably one can deter- 
mine SFRs from UV or IR alone. 

SFRs derived from dust emission are based on estimates 
of the total IR luminosity (L m = L[8 - lOOOimi]). This 
measure is accurate in the nearby Universe thanks to the 
observations of IRAS and Spitzer, both of which sampled 
the wavelength range close to the peak emission of the dust. 
Spitzer was also used to observe galaxies up to z ~ 2 in the 
mid-IR, however extrapolation to the total IR luminosity 
remains uncertain and relies on relations determined from 
nearby galaxy populations. By o bserving in the rest- frame 
far-infrared of galaxies, Herschel ( Pilbratt et al. |2010| || 1 |al- 
lows us to measure accurately this IR bolometric luminosity 
over a continuous and wide range of redshift. 
The main issue when using UV emission to estimate SFRs is 
the effect of dust attenuation. The Lir/Ltjv ratio has been 
identified as a very powerful estimator of dust attenuation 



in star forming galaxies (e.g. Gordon et al. 2000 Buat et 



alone must be used when UV and IR rest-frame data are 
not simultaneously available. Mcurcr, Hcck man, fc Calzetti| 
(1999) found a relation between the slope of the rest-frame 
UV continuum /3 (defined as f\(eig cm - 2 s _1 rrm _ ) oc A 
for A > 120 nm) and dust attenuation traced by Lir/Luv 
for local starburst galaxies observed by IUE and IRAS. This 
local starburst relation is also widely used to estimate dust 



attenuation in UV-selected galaxies at high redshift (Bur 



garella et al.||2007| |Reddy et alj[2008| . In the local Uni- 
verse luminous IR galaxies (LIRGs with Li R > 10 11 L©) are 



found to roughly follow the local starburst law ( Takeuchi et 
al. 2010 Howell et al.||2010l; normal forming galaxies, less 



active in star formation than starburst ones, do not follow 
this relation, and the spectral slope j3 is found to depend 
on a number of parameters that are not solely related to 
dust attenuation (e.g. Buat et al.|2005 Boissier et al.|200~ 



Kong et al. 20041. Checking the validity of (3 as a tracer 



of dust attenuation on large samples of galaxies at different 
redshifts is therefore a key to both understanding dust at- 
tenuation processes in galaxies, and to being able to correct 
accurately for extinction. 

In this paper we will analyse the IR and UV properties of 
galaxies over the redshift range, < z < 1. This is the 
first time accurate estimates of the total dust emission have 
been determined over a continuous range of redshift by us- 
ing several IR measurements including new Herschel data 
and which can be compared to the observed UV emission. 



al. 20051. Dust attenuation diagnostics based on UV data 



* E-mail: veronique.buat@oamp.fr 

1 Herschel is an ESA space observatory with science instruments 
provided by Principal Investigator consortia. It is open for pro- 
posals for observing time from the worldwide astronomical com- 
munity 



2 DATA 

Observations of the Lockman field were made with Herschel 
as part of the Herschel Multi-Tiered Extragalactic Survey 
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(HerMES- |p liver et al.|(|2010|))P| This field was also observed 
by GALEX jMa,rtm et al.||2005| ) in FUV (A ~153 nm) and 
NUV(A ~ 231 nm) as part of the Deep Imaging Survey 
(DIS). Only sources with z < 1 are detected by GALEX 
due to absorption by the Lyman Break. 
Our first parent sample is based on the shallow observations 
of the Lockman SWIRE fie ld (lOh 45m 00s + 58d 00m 00s, 
218' x 218' ) with SPIRE ((Griffin et al.l|2010l. We use the 



HerMES cross-identified catalogue of Roseboom et al. ( 2010 1 
which is based on a linear inversion method using the posi- 
tions of sources detected in the Spitzer 24 /jm surveys. We 
select 7435 sources, that are detected at 5<r or above at 250 
/xm, corresponding to a limit of 23 mjy (jRoseboom et aL 



20101. A redshift has been assigned to 3862 sources. Where 



available, we use spectroscopic redshifts (taken from the 
NASA/IPAC Extragalactic Database), otherwise photomet- 
ric redshifts from the SDSS/DR7 (|Abazajian et al.|2009[ ) and 
the SWIRE survey ( Rowan- Robinson et al.||2008 l are used 
by order of priority. We identify 3824 galaxies with z < 1, of 
which 30% have spectroscopic redshifts. This sample is then 
cross-matched with the GALEX detections in the NUV with 
a tolerance radius of 2 arcsec (based on Spitzer /IRAC coor- 
dinates): 2426 galaxies are detected in NUV with z < 1. The 
UV slope j3 can be calculated from the FUV — NUV colour 
for sources with z < 0.3. Of the 1542 objects at z < 0.3, 
1250 are detected in the NUV and 943 have both FUV and 
NUV fluxes that are of sufficient accuracy (Am < O.lmag) 
to enable reliable estimates of p. 

Although ~ 40% of the Lockman SWIRE sources selected 
at 250 /im are not detected in the NUV it is quite difficult 
to study the non detections since not less than 21 GALEX 
fields with different depths are needed to cover the HerMES 
field. Therefore we also consider the Herschel deep observa- 
tions of the Lockman North field ( 10h46m00s +59d01 m00s, 
35' x 35', 5<r limit at 250 /im corresponding to 8 mjy (Rose- 



|boom et aLpoTo] |): 70% of this field is covered by only 3 
GALEX fields and a check of each individual source by hand 
is therefore possible. A second advantage of the Lockman- 
North field is that it was also observed with PACS ( |Poglitch| 
et al.pOlO l and data at 110 and 160 /im are available. We 
again use the cross-matched catalogue of |Roseboom et al.| 
(2010) and perform a selection and a GALEX cross-match 
similar to that made for the Lockman SWIRE field. Photo- 



metric redshifts were determined by Strazzullo et al. (20101 



from visible, NIR and IRAC data. 246 sources are selected at 
250 /im, of which 129 galaxies are detected in the NUV. The 
58 sources that are considered as non-detections in the NUV 
arc assigned a limiting magnitude of NUV= 24.4 (Morissey 
et al.||2005 l. We consider as UV upper limits, objects with 
no GALEX source within 6 arcsec of the SPIRE position 
(for distances less than 6 arcsec the full NUV PSF falls into 
the 250/mi one)). 

Total IR luminosities, Lib., are calculated using the CIGALE 
code (|Noll et al.|2009|l wit h |Dale fc Helou| ( |2002} and |Sieben-] 
morgen fc Krugelj ( ]2007[ ) IR templates. A Baycsian anal- 
ysis is performed to deduce physical parameters from the 
fitting of spectral energy distributions of galaxies. We re- 
strict the analysis to the thermal IR (^ 24/im data), and 
defer an analysis of the full UV-to-FIR SEDs to a later pa- 
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Figure 1. Lir/Ltjv vs. Lir for sources detected in the Lock- 
man North in the redshift range < z < 1 (the redshift is 
colour-coded). Lower limits of -Lir/Luv are plotted as triangles 
for galaxies not detected in UV. The filled circles and solid line 
correspond to the relation found for local IR selected galaxies 
l |Buat et al.||2007a[ l, the empty squares and dashed line to the 
mean values per bin of IR luminosity obtained in the present 
study. Errors bars correspond to the scatter of the data (1 rms). 
The A(UV) scale is obtained using the calibration of Buat et 
|al.| ( |2005[ > . The luminosity bins and mean redshift in the bins 
are: log(Lni) < 11 and < z >= 0.4, 11 < log(Li R ) < 11.5 
and < z >= 0.55, 11.5 < log(L IR ) < 12 and < z >= 0.8, 
log(ii R ) > 12 and < z >= 0.9. 



per. For each source we consider all the available data be- 
tween 24 and 500 /j,m: by construction we will always have 
at least two data points. We keep only objects for which 
the reduced chi-squared of the fit between the best template 
and the dataset is lower than 10 (98% of the initial sam- 
ples). Typical errors on estimated Lir are found lower than 
0.1 dex. UV luminosities Luv are derived at 153 nm rest- 
frame (FUV) by interpolating FUV and NUV fluxes, a mean 
FUV — NUV colour is used when FUV data are not avail- 
able, and are defined as the quantities vL v . Both Lir and 
Luv are expressed in solar units. All magnitudes are given 
in the AB system. We assume that Q m — 0.3, S1a = 0.7, and 
H = 70 kms" 1 Mpc -1 . 



3 DUST ATTENUATION AS TRACED BY 

Lir/Luv 

The ratio of Lir/ Luv is a powerful measure of dust atten- 
uation and is known to correlate with the bolometric lumi- 



nosity of galaxies in the nearby universe (e.g. Buat et al 
|2007a|>. Spitzer observations extended the analysis out to 



0.7 (Zheng et al. 2007 Buat et al. 2007b Xu et al 



20071. Here we re-investigate this relation using the Lock- 
man North sample in the redshift range < z < 1. We 
choose this field in order to be able to discuss lower limits 
for galaxies not detected in UV. Shown in Fig. [T] is a plot of 
Lir/Luv as a function of Lir. The mean dust attenuation is 
A(UV) = 3± 1 mag for the whole sample (using the calibra- 
tion of Buat et al. ( 2005)). As expected, there is a general in- 
crease of Lir/Luv with Lir. We measure the mean Lir/Luv 
ratio within several bins of log(LiR) (using the Kaplan-Meier 
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Figure 2. Lir/-Ltjv vs. /3 (the slope of the UV continuum) for 
galaxies in the Lockman SWIRE field that are also detected in 
FUV and NUV (z < 0.3): dots represent galaxies with Ljr < 
fO 11 L0 and a robust measure of /3; filled squares represent LIRGs 
( Lih, > 10 11 Lq) with a robust measure of /3, and LIRGs with less 
reliable values are represented by empty squares. The redshift 
of sources is colour-coded. The superposed lines trace the relation 
found for local starburst galaxies (solid line) and for local normal 
star-forming galaxies (dashed line). 



estimator to account for lower limits). The locus of galaxies 
at low z appears to be consistent with the z — relations 
I Buat et al.||2007a l, however at z > 0.5 Lir/Ltjv is found 
to be lower by ~ 0.3 dex. This difference at higher z is 
of low statistical significance (of the order of the 1 a error 
bars) and could be affected by the lower-limit corrections 
but if it is real this implies a decrease in dust attenuation 
by ~ 0.7 mag in the FUV for distant, luminous IR-selected 
galaxies which is consistent with dust attenuation found for 
LIRGs at z — 0.7 by|Buat et al.|(|2007b"|). A decrease in dust 



attenuation is possibly linked to a lower metallicity in high- 
redshift systems. Indeed Reddy et al. (20101 found a corre- 



lation between oxygen abundance and dust attenuation in 
z ~ 2 star forming galaxies. The morphological evolution of 
LIRGs from mergers to more normal galaxies as z increases 



can also imply a decrease of dust attenuation (e.g. Buat et 



al. 2007b and references therein). Because of a larger gas 
mass fraction in distant galaxies, the star formation could 
be spread out over larger, more diffuse and less extinguished 
regions than in local LIRGs. A decrease of dust attenuation 
was also reported by Daddi et al. ( 2007[ ) for Ultraluminous 
IR galaxies between z = and z ~ 2 associated to a longer 
lived star formation for the higher redshift sources. 



4 DUST ATTENUATION AS TRACED BY 
THE UV SLOPE 

A common way to estimate dust attenuation in the UV in 
the absence of IR data is to use the slope /3 of the UV contin- 
uum. To estimate /3 from the FUV — NUV colour we focus 
on galaxies with z < 0.3 in the Lockman SWIRE sample and 



apply the Kong et al. ( 2004 1 recipe. A measurement error of 
0.1 mag in FUV and NUV translates into an error of 0.3 for 
j3. Therefore we restrict the analysis (otherwise specified) 



to sources measured in FUV and NUV with an accuracy 
better than 0.1 mag. We have checked that considering the 
whole sample of galaxies detected in FUV and NUV does not 
change the results, only adding a larger dispersion in the di- 
agrams. Shown in Fig. [2] is a plot of Lir/Luv as a function 



of /3, together with local relations for starbursts (Kong et al 



20041 and optically selected star-forming galaxies (Boissier 
et al.||2007| ) are superposed. Galaxies from our IR-selected 
sample exhibit a wide range of /3 and Lir/Ltjv values as al- 



ready found in the nearby universe (Buat ct al. 2005 Seibert 
et al.|2005||Takeuchi et al.|2010| >. The majority of the galax- 
ies lie between the two relations i.e. have a dust attenuation 
(as defined by Z/ir/I/uv) for a given f3 that is larger than 
that found in local normal star-forming galaxies but lower 
than in local starbursts. From the 164 LIRGs of our initial 
sample of 1542 sources at z < 0.3 in the Lockman SWIRE 
field, 125 are detected in NUV and FUV, of which 78 have 
a reliable estimate of j3 (Am < O.lmag in FUV and NUV) 
and are shown in Fig. M 67 of these 78 LIRGs (i.e. 87%) are 
found to lie below the local starburst relation. The fraction 
drops to 77% for the 125 LIRGs detected in both FUV and 
NUV without any restriction on the measurement accuracy. 
11% of the galaxies (including 10 LIRGs) detected in NUV 
have no FUV detection (for these we adopt a limiting FUV 
magnitude of 24.8 (Morissey et al. 20051 ; if we consider 



only sources with NUV magnitudes with ANUV < 0.1, all 
these sources are found below the starburst law. This esti- 
mate might be uncertain given the large number of GALEX 
fields with different exposure times. As quoted in Section 2, 
23% of the initial sample at z < 0.3th is detected neither 
in NUV nor FUV: the limiting magnitude of NUV = 24.4 
roughly corresponds to a limit on Z/ir/Ltjv of ~ 200 in the 
redshift range < z < 0.3 (for the lowest IR luminosities at 
each redshift). Therefore we miss at most 23% of galaxies 
with log(LiR/Z/uv) > 2.3 (note that fewer than 23% may be 
genuine non-detections as discussed for the Lockman North 
field) . A detailed discussion of U V non detections is deferred 
to a future paper. Nevertheless, even if we account for galax- 
ies not detected in the NUV and FUV, a substantial frac- 
tion of all the IR-selected galaxies is found below the local 
starburst law. Dust attenuation in these galaxies is overes- 
timated if the recipe for starburst galaxies based on the UV 
spectral slope is applied. 

The most distant galaxies of our sample, including LIRGs, 
are found to depart strongly from the starburst law (see 
Fig. [2}. j3 is calculated from the observed FUV - NUV 
colour: the wavelengths sampled go from 150 and 230 nm 
at z=0 to 118 and 180 nm at z=0.3. Therefore, the power 
law model for the UV continuum assumed to calculate p 
might not be valid for the whole 118-230 nm range and for 
dust attenuations ranging from 0.5 to 5 mag in UV. The un- 
certainty can also be amplified because of the rather large 



bandwidths of the GALEX filters (Morissey et al. 2005) 



Variations in dust properties and star versus dust geometry 
are known to induce large dispersions in the -Lir/Ltjv — 
diagram ( Gordon et al.|2000 Boquien et al.|2009 |). Our se- 
lection at 250 fim may also lead to more quiescent galaxies 
than local LIRGs in terms of star formation which is known 
to increase /3 for a given Lir/Luv (Kong et al. 2004). At 
z=2 |Reddy et al.| ( |2010[ > found that Lyman Break Galaxies 
also detected at 24 /jm and with bolometric luminosity lower 
than 1O 12 L0 follow the local starburst law, although with a 
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Figure 3. Fraction of SFR measured in IR, SFRjr as a function 
of the IR luminosity, Ljr. The mean and standard deviation of 
the values in 5 bins of log(L IR : < 10, 10-10.5, 10.5-11, 11-11.5, 
> 11.5 are shown in red, joined by a solid line to guide the eye. The 
blue dashed line represents the fraction of galaxies detected in the 
UV at each value of log(Ljjj) 



large dispersion. These galaxies have f3 < —0.5 and therefore 
are much bluer than the bulk of our sample. At high red- 
shift the dynamical range in stellar population is narrower 
and the contribution to the UV continuum of stars not re- 
lated to the current star formation is likely to be negligible 
in star forming galaxies at z > 1. The spatial distribution of 
dust and young stars may also be similar to that found in 
strong local starbursts making the law proposed by |Meurer,| 
Heckman, & Calzetti ( 1999[ ) more appropriate for these high 
redshift objects than for our present sample. 

This topic will be re-investigated in a future work by 
fitting the whole SED from the UV to the IR both at low 
and high redshift. 



5 CALIBRATING MEASURES OF STAR 
FORMATION RATE 

An important check is to test whether IR emission alone 
provides a robust estimate of the total SFR over the whole 
range of IR luminosities or if ignoring the direct UV light 
leads to systematic errors. We again use the sample derived 
from the Lockman-SWIRE field because of the large number 
of galaxies in the sample. 

We take SFRir, + SFRtjv = SFR to t as a reference for the 
SFR measure where Iog(SFRm) M yr -i = log(Z/iFt)L —9.97 
and log(SFRuv)M yr- 1 = log(Luv)L —9.69. The UV emis- 
sion is not corrected for dust attenuation and we use the SFR 



calibrations from Buat et al. ( 2008 1 with the assumption of 
a constant SFR over 10° years and a Kroupa initial mass 
function ( Kroupa|2001 1. 

Plotted in Fig. |3j is the SFR fraction measured by the IR 
luminosity, SFRiR./SFR to t, as a function of Lib,, with an 
additional curve showing the fraction of the galaxies in the 
sample with that IR luminosity that are detected in the 
UV. When the sample is considered as a whole, we find 
that SFRir, measures ~ 90% of the total SFR. This fraction 



varies with Lm,: from 94±10% for the most luminous galax- 
ies of our sample (Lib. > 10 11 ' 5 L Q ) with 45% detected in 
UV, down to 71 ± 17% when Lir < 10 10 Lq or equivalently 
SFRtot < IMoyr -1 and a UV detection rate reaching 85%. 
This demonstrates that the combination of obscured and 
unobscured SFRs is required to determine accurate SFRs in 
galaxies with low star formation activity. Intrinsically faint 
galaxies detected in rest-frame UV surveys are found to be 
a significant component of luminosity functions and of the 
total star formation density at both low and high redshifts 



( |Buat et aL 2007a Reddy fc Steidel|2009[ ). As a consequence, 
both UV and IR selected samples must be built and their 
contribution added to measure the total star formation den- 



sity at a given redshift ( Iglesias-Paramo et al. 2006 Reddy 
fc Steidel||2009l 



Shown in Fig 111 is a comparison of the SFR measured in the 
UV with SFRtot, for both uncorrected and corrected UV 
luminosities, where the correction used is that of |Meurer^1 
Heckman, & Calzettil (19991): A uv = 4.43 + 1.99,0 deduced 



from their Ltr/Ltjv-/3 relation. For this correction we only 
consider galaxies with z < 0.2 to narrow the wavelength 
range used to calculate /3. Without dust attenuation correc- 
tions, SFRtjv underestimates the total SFR by a factor ~ 6 
for SFRtot ^ 1 Mayr" 1 and ~ 30 for SFR tot ~ 100 Meyr" 1 . 
As expected from Fig [2] with dust corrections based on (3, 
SFRs are ov erestimated by a factor ~ 2 — 3. At z ~ 1 Elbaz 
et al.|(|2007|) only obtain SFR s; 10 M yr _1 when correcting 



UV luminosities for dust attenuation based on values of j3. 
They use the U — B colour to measure f3 which corresponds 
to the wavelength range 180-225 nm at Z ~ 1; which is dif- 
ferent from the ranges used in this work and by |Meurer,| 
Heckman, fc Calzetti| (|1999[). |Burgarella et al.| (|2007[) and 



Reddy & Steidel ( 2009 I corrected UV measurements using 



the Lir/Ltjv-/3 relation for local starbursts and found SFR 
deduced from UV-corrected luminosities to be in agreement 
with SFRtot out to few 100s M Q yr _1 , allbeit with signif- 
icant dispersion. As discussed in section 4, this agreement 
is likely to be due to the different nature of UV selected 
galaxies at z > 1 as compared to the IR selected galaxies 
analysed in the present work. 
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Figure 4. A plot of SFR UV vs. SFR tot , where SFR tot = SFR IR + 
SFRtjv- SFRs deduced from uncorrected UV fluxes are plotted 
with blue dots, the largest ones correspond to galaxies at z < 0.2. 
SFRs deduced from UV fluxes corrected for dust attenuation with 
the relation of |Meurer, Heckman, &: Calzetti| ( |1999fr are plotted 
with open green circles for galaxies with z < 0.2. The linear re- 
gression lines between SFRtot and the two estimates of SFRtjv 
are plotted as dashed lines. The black solid line corresponds to 
equal values on both axes. 
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